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Введение. Важные детали машин и механизмов, испытывающие значительные нагрузки, изго- 
тавливаются из сталей, и их поверхность упрочняется такими технологическими операциями, как 
поверхностная закалка токами высокой частоты, химико-термическое упрочнение, плазменно- 
диффузионное, лазерная обработка, упрочнение наклепом, нанообработка и др. Для обеспечения 
высококачественного производства необходимо иметь различные средства неразрушающего 
(в некоторых отраслях — стопроцентного) контроля, которые позволят не только выявлять и изо- 
лировать брак, но и проводить диагностику причин появления брака для корректировки техноло- 
гического процесса. 

Применяемые методы настройки приборов для неразрушающего контроля, использующие 
электромагнитные датчики, связаны с изготовлением и экспериментальным исследованием на- 
строечных образцов с известными (а зачастую и заданными) свойствами. Причем количество этих 
образцов должно быть статистически значимым (иногда более 100 штук). Изготовление такого 
количества образцов и исследование их физико-механических свойств требует как больших фи- 
нансовых затрат, так и затрат времени, что делает невозможной быструю переналадку системы 
контроля на новую задачу. 

Проблему можно решить, если применить компьютерное моделирование [1]: математиче- 
ская модель описывает связь между сигналом датчика и физико-механическими свойствами изде- 
лия, а имитационная модель позволяет создавать «виртуальные» образцы изделия с требуемыми 
свойствами. Таким образом, методом компьютерного эксперимента возможно определять опти- 
мальные условия контроля и находить значения коэффициентов для настройки приборов нераз- 
рушающего контроля. 

Физико-математическая модель образования сигнала датчика. В приборах неразрушаю- 
щего контроля в качестве накладного электромагнитного преобразователя часто используется 
датчик, показанный на рис. 1. 

Изделие с переменными по толщине (координате Х’) свойствами представляется состоя- 

щим из и тонких слоев, имеющих неизменные электромагнитные свойства: толщину й, комплекс- 


ную магнитную проницаемость п, электропроводность о, коэффициент Рэлея у. Датчик состоит 


из ферромагнитного сердечника с намотанными на нем возбуждающей обмоткой И’, и измери- 
тельной обмоткой И/>. Он расположен над контролируемым изделием с зазором 6 [2]. 


Электромагнитное поле в проводящей среде описывается уравнениями Максвелла: 
гоёН =оЕ г гоЕЁ = —тоиН Е 
Для схемы расположения датчика над контролируемой средой можно записать граничные 
условия: 
ОН Е 
—=0, —=0, Н,=0. 
бу 7 
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Рис.1. Электромагнитный датчик, расположенный 
над многослойной контролируемой средой 


Следовательно, 
| ан. дн ЗЕ 
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Получаем уравнение: 
: + ге — юноЕ. =0, 
ду ОХ 
которое является однородным уравнением Гельмгольца в декартовой системе координат: 
0°Е. 0?Е. 
дх* у ду” 








+т?Е. =0, 








где т” = —]оно. 
Дальнейшее решение уравнения позволяет определить ЭДС датчика [3]: 


|1 154 С7” Ат 400’. ‚м 
Е =-Й И’, [о И ЧЕ ТЕ о УИ 
йа 2 ВАЙ Пе ВА ел 
Полученное уравнение можно использовать для определения сигнала электромагнитного 
датчика в имитационной модели. Для примера приведено сравнение моделированных И измерен- 
ных на образцах изделий сигналов датчика, из которого хорошо ВИДНО практически полное их 


совпадение (рис.2). 


т] 
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Рис.2. Сравнение амплитудно-частотных характеристик моделированного и измеренного 
на образцах изделий сигналов при различных относительных зазорах между датчиком и изделием 
В ряде случаев, применение основных колебаний для электромагнитного контроля приво- 
дит к недопустимой погрешности. Например, ЭДС датчика в зависимости от толщины упрочненно- 
го слоя определяется, главным образом, отражением электромагнитной волны от поверхности, 
разделяющей упрочненный слой и основу. Упрощенно зависимость ЭДС от толщины упрочненного 
слоя в этом случае определяется коэффициентом отражения: 


-1, (*) 





где и! и с! соответствуют упрочненному слою, а и> и о> — основе. 
При всех известных видах упрочнения, вызывающих физические дефекты в идеальном 
кристалле железа, магнитная проницаемость д, и удельная электропроводность с, упрочненно- 


го слоя уменьшаются по сравнению с магнитной проницаемостью п, и удельной электропровод- 
ностью о, основной стали. Если и и с изменяются пропорционально, то выражение (*) стре- 


мится к нулю, и преобразователь теряет чувствительность к толщине упрочненного слоя. В дей- 
ствительности выражение (*) не равно нулю, но близко к нему. Сигнал отраженной волны на 
границе раздела слоев настолько мал, что его превосходят помехи, вызванные вариациями элек- 
тромагнитных свойств изделия. Это и является причиной возникновения недопустимо большой 
ошибки при контроле некоторых видов упрочнения по основным гармоникам сигнала преобразо- 
вателя. 

Исследование зависимости информативности комбинационных колебаний по сравнению с 
основными показало, что в ряде случаев комбинационные колебания значительно информатив- 
нее. Несмотря на то, что на комбинационные колебания сильно влияет зазор между датчиком и 
контролируемым изделием, влияние толщины упрочненного слоя на них гораздо сильнее. Поэто- 
му отношение полезного сигнала к помехе в них больше, чем в основных гармониках. Использо- 
вание комбинационных колебаний повышает точность контроля толщины упрочненного слоя. 

Например, при упрочнении борированием коэффициент отражения [см. формулу (*)] ра- 
вен 0,058, тогда как аналогичный коэффициент для комбинационных колебаний равен 0,44 - это 
почти в десять раз выше. 


1059 


Технические науки 


ЭДС комбинационных колебаний в слабых полях рассчитывают, учитывая коэффициент 
Рэлея х и то, что электромагнитные свойства слоев изделия изменяются по глубине (см. рис.1). 
Решение двумерного дифференциального уравнения Гельмгольца для каждого 1-го слоя 
модели имеет вид: 
н„=— УИ, екр-Лх) — Си окр ров Ху. 
ОН; 
Граничные условия между слоями 1 и 1+1 можно представить выражениями: 
Г: 2 © и] 
= | Авехр(- ЛА.) — Си ехрО 9х, |= 
И; 
Е | .) (х) 5) (>) 
р | Ал ехр(-— 7 +1АХ,) ет. ехр( А их, | ) 
И 
Ау ехр(-—//А92х,) те ее ехр( Арх, = 


Е Аня ехр(-— ЛА их, р фт ехр( АР х,). 


В последней (и+1)-й области металла (т.е. на бесконечно большой глубине) электромаг- 
>0иН 





нитное поле полностью затухает, поэтому при х-—*® Н 0. Такое возможно 


пу п-+1,х 


при А,1«=0. 


В модели выбирается большое количество настолько тонких слоев, что напряженность 
магнитного поля внутри каждого слоя изменяется менее чем на 1%. Большим количеством слоев 
достаточно точно дискретно отражается непрерывное распределение электромагнитных СВОЙСТВ 
упрочненного слоя и основного металла. 

Напряженность поля в указанных точках определяется как функция времени: 

Ни =Ни (Е -+уи;). 


Приведенная последовательность повторяется дважды для двух частот © И ©, основных 


колебаний тока возбуждения преобразователя, которые в нелинейной среде, взаимодействуя, 
образуют комбинационные колебания. 

Чтобы комбинационные колебания отличались по частоте от основных колебаний и их 
высших гармоник, основные частоты должны относиться как достаточно большие целые взаимно 
простые числа № и М: 

61:0. =ЛМ: М. 

Напряженность результирующего поля в точке х,: 

Н;= Нид: (ФИ +) + Н и, (Е + у, ). 

Получившееся колебание -— периодическое; его период составляет 

2п 21 
Г=—М=— М. 
0 (> 

Непрерывно изменяющаяся во времени напряженность вычисляется дискретно через ин- 
тервалы времени ЛЕ. Если ограничить спектр сигнала суммарными комбинационными колеба- 
ниями третьего порядка, то, по теореме Котельникова, интервал времени 

Г 


р: 
2М+АМ 
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Индукция магнитного поля в точке х, в момент времени {, определяется уравнением 
Рэлея: 
В =НН» + х.НьяепН,, | 
Рабочий поток магнитной индукции в контролируемом изделии рассчитывается заменой 
интеграла суммой: 


п-+1 


Ф, = В => Врй,, 
0 1=1 


где Ф, есть дискретная функция времени. 


Гармоническим анализом переменного магнитного потока можно определить веществен- 
ную и мнимую составляющие потока суммарных комбинационных колебаний третьего порядка: 


Т т 
а | Ф(В (2, +) = УФ, о А 
78 ' т‘ т 


ком 


Т 
9 2х 2л 
Ф/ = = | Ф@соз(2, +01) =— > Ф, с05——(2М + М№р, 
0 р 


где т- число интервалов времени ЛЕ, на которые разбит период. 
Комплексная амплитуда рабочего потока: 
Ф ом — Ф ом о 
частично сцеплена со вторичной обмоткой преобразователя, поэтому комплексная ампли- 
туда комбинационных суммарных колебаний третьего порядка рассчитывается следующим 
образом [4]: 
Эм =-Л®, +2®,)\.,Ф, . 

Подобным образом с помощью гармонического анализа магнитного потока в математиче- 

ской модели вычисляются ЭДС основных колебаний и разностных комбинационных колебаний 
третьего порядка. Принципиально возможно и вычисление высших гармоник. 
Имитационная модель создания виртуальных образцов. Предназначена для имитации 
виртуального прибора для неразрушающего контроля, виртуальных образцов, ЭДС сигналов из- 
мерительного преобразователя (датчика), вычисления линейной регрессионной модели по сигна- 
лам основных и комбинационных колебаний, а также определения ошибки контроля на виртуаль- 
ных проверочных образцах. 

Имитационная модель реализована в виде программы для ЭВМ. Структурная схема имита- 
ционной модели показана на рис.3. При моделировании вычисляются ЭДС датчика для набора 
значений либо зазора, либо толщины упрочненного слоя. Кроме того, при вычислении ЭДС 
происходит случайное варьирование электромагнитных параметров виртуального образца изде- 
лия, а также случайное варьирование того параметра, который не варьируется детерминирован- 
но. В конце моделирования к полученным значениям ЭДС добавляются случайные значения, ими- 
тирующие аддитивную и мультипликативную погрешности прибора для неразрушающего контро- 
ля. Гакие вычисления ЭДС датчика производятся для заданного ряда частот как основных, так и 
комбинационных колебаний. 

На рис.4 показана структурная схема блока вычисления ЭДС сигналов измерительного 
преобразователя при детерминированном варьировании зазора между датчиком и изделием. Ана- 
логичный вычислительный блок используется при детерминированном варьировании толщины 
упрочненного слоя, тогда производится цикл не по значениям зазоров, а по значениям толщин 
упрочненного слоя. 
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Рис.3. Структурная схема имитационной модели 


В системе создается два набора виртуальных образцов — настроечные и проверочные. На- 
строечные образцы используются для построения регрессионной модели. После получения коэф- 
фициентов уравнения регрессии, в это уравнение подставляют значения сигналов, полученные на 
проверочных образцах. То есть в этом случае имитируется процесс измерения контролируемого 
параметра виртуальным прибором для неразрушающего контроля. Результат этой операции - 
показания прибора, получающиеся при измерении образца, не участвовавшего в настройке. Оба 
набора виртуальных образцов могут различаться между собой. 
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Рис. 4. Структурная схема вычислительного блока ЭДС 
сигналов измерительного преобразователя 


При моделировании возможно задавать различные начальные условия для массивов на- 
строечных и проверочных образцов. Например, количество настроечных образцов можно зада- 
вать небольшим, чтобы оценить представительность выборки настроечных образцов в данных 
условиях контроля. Различной может быть и величина случайных вариаций электромагнитных 
свойств в каждом из наборов. Тогда появляется возможность определить, как повлияет ошибка в 
оценке параметров настроечных образцов, или что покажет прибор при контроле изделий из дру- 
гой партии, отличающихся либо значениями технологических факторов, либо их разбросом. 

Ошибка контроля зависит от количества настроечных образцов при различном количестве 
частот, составляющих спектр измерительных сигналов. С ростом количества настроечных образ- 
цов ошибка сначала падает, а потом перестает меняться. Кроме того, ошибка контроля зависит и 
от числа частот, составляющих спектр сигнала датчика. При малом количестве настроечных о6б- 
разцов с ростом числа частот, составляющих спектр, включенный в регрессионную модель, растет 
и ошибка контроля. Это происходит из-за снижения количества степеней свободы при построении 
модели и, следовательно, снижения представительности выборки сигналов. 
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Имитируемый образец разбивается на и слоев с неизменными электромагнитными свойст- 
вами внутри слоя. Магнитный поток в изделии вычисляется интегрированием по толщине изде- 
лия. Поэтому толщина одного слоя разбиения должна быть достаточно мала, чтобы изменение 
величины напряженности поля при переходе от слоя к слою также было малым. Однако при 
слишком большом числе слоев резко растет вычислительная ошибка и продолжительность вы- 
числений, следовательно, требуется оптимизация числа слоев для минимизации погрешности 
моделирования. 

Ошибка контроля зависит и от числа слоев, на которые разбивается виртуальный образец. 
При количестве разбиений более 1000, ошибка контроля начинает расти. Это происходит из-за 
накопления погрешностей округления при вычислениях с конечной точностью [5]. 
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Рис.5. Частотные характеристики модуля передаточного коэффициента 
толщины упрочненного слоя для комбинационных колебаний 


Для обеспечения высокой точности измерения контролируемого параметра должен быть 
установлен оптимальный спектр частот возбуждения преобразователя. Амплитудно-частотные 
характеристики передаточного коэффициента толщины упрочненного слоя (рис.5) позволили 
установить области частот возбуждения датчика, при которых он высокочувствителен к контро- 
лируемому параметру, и выбрать, таким образом, составляющие многомерного сигнала. 
Заключение. Разработана математическая модель образования сигнала электромагнитного дат- 
чика, расположенного над многослойной проводящей ферромагнитной средой для основных и 
комбинационных колебаний. Несмотря на то, что амплитуды комбинационных колебаний намного 
меньше, чем основных, они оказываются более информативными. Используя комбинационные 
колебания, можно существенно повысить точность неразрушающего контроля. Поскольку при 
упрочнении значительно изменяется коэффициент Рэлея, то комбинационные колебания оказы- 
ваются намного чувствительнее к толщине упрочненного слоя, чем основные. 

Разработана имитационная модель образования сигнала электромагнитного датчика для 
основных и комбинационных колебаний. Поскольку аналитическое выражение для комбинацион- 
ных колебаний не очень удобно для прямого расчета, то предложен рекуррентный метод, исполь- 
зуя который можно достаточно просто получить значения ЭДС датчика как для основных, так и 
для комбинационных колебаний (а также для высших гармоник, если такое понадобится) любого 
порядка. Генерация виртуальных образцов и их использование для построения и проверки точно- 
сти регрессионной модели позволяет определить оптимальные условия контроля без изготовле- 
ния и исследования реальных образцов изделий. Кроме того, значительно сокращается количест- 
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во реальных образцов изделий, требуемое для окончательной настройки прибора при его экс- 
плуатации. 

Усовершенствованы методы формирования и обработки сигнала датчиков в целях повы- 
шения точности контроля. Использование имитационной модели позволяет определить оптималь- 
ный спектр возбуждения датчика и параметры регрессионной модели для обработки его сигнала. 
Выбор оптимального спектра возбуждения значительно повышает точность контроля и уменьшает 
влияние на его результаты случайных вариаций свойств изделия, а также различных помех, воз- 
действующих на устройство в целом. 
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